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1.0 STATIC STRATIFICATION FOR SIGNATURE EXTENSION 


1.) OBJECTIVE 

The objective of this task is to evaluate and If necessary Improve 
the signature extension stratification maps developed by UCB In previous 
LAC IE tasks. Specifically, the ability of the strata to group spectrally 
similar wheat subclasses Is being evaluated. In addition, variables used 
to produce the strata maps are being analyzed with regard to their Influence 
on spectral signature of wheat fields, and, If needed, new variables will 
be introduced to refine or otherwise alter the stratification process, 

1.2 GENERAL APPROACH 

The task, designed to evaluate the statistical significance of the 
static stratif fcatlon and to provide Information for Its refinement, has 
been divided Into two subtasks. These are (1) to evaluate through 
analysis of variance (ANOVA) the static stratification’s ability to group 
spectrally similar areas In order to maximize signature extension success 
and (2) to determine the statistically significant signature controlling 
variables for use In refining the stratification procedure. 

SUBTASK A; SPECTRAL GROUPING ANALYSIS 

The purpose of this subtask Is to determine If the static stratification 
does In fact Isolate areas tending to have similar wheat signatures. This 
analysis Is also Intended to discover the extent to which Individual strata 
can be grouped together and still provide for potentially successful 
signature extension. 

The Initial experimental procedure consisted of obtaining a systematic 
sample of wheat fields In each land use/sol Is/cl Imatic stratum and then 
comparing, thi ough analysis of variance, the differences In spectral 
signature for the corresponding fields within and between strata. Results 
from this analysis (reported In Thomas and Hay, 1977) suggested the need 
for a revised approach. Rather than compare differences In signature 
between Individual wheat fields, a procedure that compared differences (1) In 
wheat spectral density functioning between strata and (2) between wheat 
spectral density functions and those for other crops appeared more appropriate. 

This second approach Is Intended to more effectively Identify potential 
classification accuracies, and therefore ensure signature extension success, 
within and between strata. 

It is composed of four basic steps. The first, preprocessing , Is 
Intended to standardize segments to a common sun elevation and haze condition. 
This Is accomplished by Implementation of XSTAR haze correction procedures 
(Lambeck 1977) developed at ERIM, Preprocessing in this case is Intended 
to provide a more stable measurement frame (Landsat or Tassel led Cap Space) 
and thereby Increase the ease with which real spectral differences can be 
identified and evaluated. 








The second step Involves clustering via ISOCLAS (adapted from JSC) 
each preprocessed sample segment. ISOCLAS control parameters are patterned 
after those currently In use at JSC, adjusted when necessary to optimize 
clustering performance based other UCB LACIE experience. Thirdly, each 
resulting cluster, or combination of similar clusters, Is then label led as 
to cover type. This Is accomplished by drawing a random sample of pixels 
for each segment to a pre-specif led sample size fy; each cluster. An 
Image analyst then references each pixel drawn to the correspond Ing ground 
data (JSC Blue-ray sheets) and labels the pixel as to cover type. When 
ground data Is missing or questionable, the analyst Interprets multidate 
PFC transparencies to Identify the cover type. 

Data for all labelled pixels In all clusters Is summarized by land 
use/sol 1 assoclatlon/cl Imatic stratum In each sample segment. The 
proportion of each cover type (e.g. wheat, alfalfa, corn, fallow, etc.) 

In each cluster Is also determined from the pixel sample. Clusters are 
then assigned to cover type groups according to the dominant cover type 
In each. The range In spectral reflectance mean vectors (l,e. vectors 
containing the average response In each handsat or Tassel led Cap band) 
for clusters belonging to a given cover type group Is defined to represent 
the range of the spectral reflectance density function corresponding to 
that group. Multiple comparison tests (based on ScheffI 1959) are then 
applied to determine the statistical separability of each cluster belonging 
to wheat cover types relative to non-wheat clusters In the same as well 
as different signature extension strata. This comparison utilizes the 
spectral mean vectors and covariance matrices for the cluscers as defined 
In the clustering step. Results are Interpreted in terms of the degree of 
overlap between cover type spectral density functions. Strata are combined 
when their respective wheat spectral density functions have significant 
overlap and overlap with other cover types Is minimal or Is similar In 
both strata. In the limiting case, the stratification should give rise 
to distinctly different wheat spectral density functions between strata. 

Two biophase periods have been Initially selected In Kansas and North 
Dakota In which to apply the revised grouping analysis procedure Just^ 
described. Date ] In both states represents a wheat einergence condition. 
The second date corresponds approximately to a jointing or advanced jointing 
condition for the wheat crop. These time periods were selected based on 
sensitivity analysis results reported In section 1.3 which suggested these 
biophase stages were most difficult to characterize by static stratification 
variables. This analysis is therefore considered conservative relative to 
the performance of the static stratification. Tables 1.1 and 1.2 list the 
sample segments employed for each date In each state. Available segments 
were limited to those having ground data to minimize Incorrect I nterpretatlon 
of results. 
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Table 1.1: Kansas San<i le SegiTicnts Currently In Use for Spectra) Grouping 

Analysis 


Seg. ID 

Date § 1 

Date tf 

1019 

1 - 19-76 


1020 


5-7-76 

1035 


5 - 6-76 

10t»l 


5-6-76 

1154 

1-18-76 


1171 


5 - 4-76 

1176 



1851 

1-19-76 


1852 


5-6-76 

1854 



1855 

1-14-76 

5 - 6-76 

1857 

1 - 20-76 


i 860 


5 - 6-76 

1861 

1 - 20-76 

5 - 7-76 

1864 

1 - 20-76 


1880 

1 - 18 , 19-76 

5 - 6-76 

1882 

1-18=76 


1884 


5-4-76 

1886 


5 - 6-76 

1887 

1 - 18-76 

5-6-76 

1891 

1-18-76 



TOTALS 


21 


12 (11 dif. segs.) 12 




Table 1.2: North Dakota Sample Segments Currently In Use for Spectral 

Grouping Analysis 


Seg. ID 

Date ^ 1 

Date ff 2 

1603 

5-28-76 

7-4-76 

1614 


7-1-76 

1618 

5-24,26-76 

6 - 30-76 

1624 

5-25-76 

1633 

5-26-76 

6 - 30-76 

1642 

5-24,25-76 

6 - 30-76 

1645 

5-24,25-76 

6 - 29-76 

1647 


7-4-76 

1648 

5-28-76 


1650 

5-27,28-76 


1651 

5 - 28 , 29-76 


1656 


7-2-76 

1660 

5 - 26-76 • 


1662 

5-24-76 



TOTALS 


14 


16 (|1 dif, segs.) 8 
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SUBTASK B; SIGNATURE CONTROLLING FACTOR SENSITIVITY ANALYSIS 

This subtask Jtffers from the first In that the cause for signature 
variability Is explored. The basic approach Is to develop regression 
relationships relating spectral reflectance (dependent variable) to a set 
of static stratification, seasonal, and date-specific (predictor) variables, 
A list of these variables Is given In Table 1.3. The Importance of the 
static (long term, I.e. 30 year plus), seasonal, and date-specific variables 
Is then expressed In terms of the percent of total Landsat band variance 
accounted for by each. Reflectance versus predictor variable regresslor^ 
are established by blophase to Identify the changing Importance of the 
variables over time. 


Two different methods of obtaining mi'tched spectral response and 
predictor variable data have been developed In this study. The first, 
completed this past quarter, was applied to a subset of Kansas 1973*’7^ 

T5E segments (see Table 1,^). Briefly, the procedure consisted of first 
sun angle-correcting Landsat data for all segments via JSC tables (Lockheed 
1976 ), generating a band 7/5 ratio for each pixel, and generating Landsat 1 
Tassel led Cap spectral representation for each segment, A small (<_ 15) 
systematic sample of wheat fields was drawn from each land use/sol Is 
stratum In each segment examined. Field-specific Landsat and Tasselled 
Cap band averages were then obtained for each field sampled using an 
Interactive color display system. Values for static, seasonal, and date- 
specific variables were measured according to procedures listed In Table 
1.3 and recorded by field along with the Landsat spectral data. 

Once all data had been obtained and coded for each field sampled, 
regression relationships were established betv^een each spectral band and 
list (Table I. 3 ) of signature predictor variables. Separate regression 
equations wc!*e obtained from Landsat bands 4, 5, and 7 and Tasselled Cap 
bands I, 2, 3, and h by date, Landsat band 6 was not Included due to 
noise problems and data storage/manipulatlon limitations present at the time. 

The relative Importance of each signature predictor variable listed 
In Table 1,3 Is expressed two ways. Measure I consisted of the percent 
of total spectral variance (by band) explained by the addition of a given 
predictor variable to the regression equation. Variables were added In 
the same order as listed in Table 1.3 using s stepwise regression technique. 
The order - static, seasonal, date-specific - was chosen to most effectively 
Identify the percent spectral variance accounted for by the static 
stratification variables before application of a signature extension 
algorithm. The R^ (multiple correlation coefficient squared) Increments, 
representing the percu.nc of variance added by variable, are highly 
dependent on this ordering. 


The second measure of signature predictor variable Importance did 
not depend on order of entry Into the regre.nslon. A forward selection 
regression procedure was used to order variables and tabulate the R^ 
Increments. Using this technique, the predictor variable having the 
highest simple correlation with two spectral band In question Is entered 
Into the regression first. The next variable entered Is the one having 
the highest partial correlation with the spectral band after the effect 
of the first variable entered Is removed. The third variable entered has 
the next highest partial correlation with the spectral response variable 
among all remaining predictor variables with the effects of the first two 
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Table 1,3: Signature Predictor Variables Used In the Kansas 1973”?^ Wheat 

Spectral Sensitivity Analysis 


Pred i ctor Variable 

Static Stratification Variables 
(Obtained from Static Strata map) 

A. Cultivated area percent 
(CULTPCT) 

B. Crop diversity 

(CDIVSITY) 

C. Soil aval lable water 
holding capacity 

(AWC) 

D. Long term average growing 
season degree-days 

(LTGSDD) 


E. Long term average growing 
season precipitation 

(LTGSP) 


Measurement Technique 
Used for Each Field Sampled 


Midpoint of cultivated area 
percent range for the land use 
class covering the wheat field. 

Midpoint of the crop diversity 
class covering the wheat field. 
Crop diversity Is coded in 
terms of three equal -si zed 
classes together spanning a 
scale of zero (lowest) to 
ten (highest). 

Average Inches of water held 
per inch of soil at field 
capacity In the top 24 inches 
for the static strata soil 
association covering the wheat 
field. These values are obtained 
from information /vailable in 
county soil survey publications. 

Midpoint of growing season 
degree-day class covering the 
wheat field. Degree-day 
classes obtained from 30 year 
average data by automatic and 
manual interpolation of ground 
meteorological station data 
for the period April through 
June . 

Midpoint of growing season 
precipitation class covering 
the wheat field. Precipitation 
classes obtained from 30 year 
average data by automatic and 
manual interpretation of ground 
meteorological data for the 
period April through June. 
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11, Seasonal Variables 

(Specific to 1973~7^ Growing Season) 

« 

A. Grovel ng season degree-days 

accumulated to Landsat pass date 

(SUMGSDD) 


B. Growing season precipitation 
accumulated to Landsat pass 
date 

(SUM6SP) 

C. Average January 197^ 
temperature 

(JANTEMP) 

D. Planting season (1973) degree- 
days accumulated to Landsat 
pass date 

(SUMPSDD) 

E. Planting season (1973) 
precipitation accumulated 
to Landsat pass date 

(SUMPSP) 

III. Landsat Date-Specific Variables 

A. Average bare soil albedo 

(ALBDRY) 


B. Precipitation In the four 
days preceding Landsat pass 
date 

(PPT4DA) 

C. 100 X Tangent of Landsat 
scan angle 

(SCANANG) 
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Calculated from temperature 
dates supplied from nearest 
ground meteorological station 
hav I ng a phy s I ca I /c 1 1 ma 1 1 c 
setting most closely approximating 
the segment In which the wheat 
field falls. Growing season 
period: April through June. 

Determined as In 11. A. 
relative to precipitation data. 


Determined from nearest 
meteorological station as 
In II. A. 



Determined as In 11. A. but for 
the period September through 
November . 


Determined as In 11. A. relative 
to precipitation data In the 
period August through November. 


Average sum of Landsat band k, 5, 
and 7 digital reflectance values 
for a sample of bare soil fields 
falling In the land use/soMs 
stratum In which the wheat field 
falls. 

Determined as In 11. A. relative 
to precipitation data. 


Departure measured along scan line 
of segment relative to an Imaginary 
base line perpendicular to the scan 
direction and pass through the 
Landsat full frame center point. 
Measurement made on full frame 
containing segment for appropriate 
pass date. The departure, reported 
In nautical miles, is defined as 
zero on ^^e base line and Increases 


E 

if 









r 




0, Average Slope* 
(SLOPE) 


E. Landsat Band 7 to Band 5 ratio 
(RASP) 
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positively to the east and 
negatively to the west. Then 


tan (scan angle) 


departure (n.m.) 
mean sa~t , a 1 1 i tude 
[kBk n.m.). 


Use 1:250,000 topography sheets 
to calculate average slope for 
the land use/soils stratum 
covering the wheat field in 
question for the segment In 
question. 

Obtain average (2x) band 7 to 
(lx) band 5 ratio for the wheat 
field. 


Footnote: Slope was fixed, not date-specific. Its placement in this 

position of predictor variable listing was for the purpose of 
determining the R^ overlay between scan angle and slope. 


I 
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. 4s Table 1,45 Number of Fields Sampled by Segment by Date for the Subset 

of Kansas 1973-74 Segments Used In the Hrst Spectral 
Sensitivity Analysis 


Date Set 


Segment 

No, 1 
One of: 

201 0 , 

Oct, 
24 J 

No. 2 
One of; 

No. 3 
One of: 

No. 4 
One of: 

26 1 May 

I7J 

. No. 5 
14 June 
‘74 

1018 


17 


17 


1026 


10 


10 


1029 

14 

14 


14 


1036 

37 


37 

37 

37 


(Land sat 
only) 




(Landsat 

only) 

1045 North 




22 


1065 

10 


10 

10 

10 


(Landsat 

only) 




(Landsat 

only) 

1109 

15 

15 

15 



TOTALS 

76 

56 

62 

1 10 

47 


(47 Landsat 
only) 




(Landsat 

only) 
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variables entered removed, and so on. Avera<je rank (order) and R Increment 
among all bands for a giben date provide the second measure of performance. 

A second method of obtaining the matched spectral response and predictor 
varlcble data will be utilized with 1976 Kansas and North Dakota LACIE segments. 
This procedure wi 1 1 differ from the first In the following respects. Spectral 
data will be obtained from a systematic or random sample of pixels labelled as 
wheat In the spectral grouping analysis described In Subtask A. This Implies, 
also, all preprocessing and clustering steps discussed in that analysis. A 
more refined set of static, seasonal, and date-specific variables will be 
evaluated for signature response correlation. These signature controlling 
variables will be selected based on the results of the first analysis. 

1.3 RESULTS LAST QUARTER 

SUBTASK A; SPECTRAL GROUPING ANALYSIS 

Work In progress. 

SUBTASK B: SIGNATURE CONTROLLING FACTOR SENSITIVITY ANALYSIS 

Ordered Regression 


Results from analysis of the subset (Table 1.4) of the 1973“74 Kansas 
Test S Evaluation Lacie segments are presented in this section. Discussion 
will focus on results for the first measure of signature predictor. variable 
importance, namely .percent of total spectral variance (by band) explained 
by addition of a given predictor variable to the regression equation.' The 
reason for this emphasis is that the ordered entry of variables - static, 
then seasonal, and then date-specific - should most efficiently identify 
the potential signature variance control to be gained with static and 
seasonal stratification prior to date-specific signature extension. 

Tables 1.5 and 1.6 summarize the Incremental proportion of wheat 
spectra^ variance (equivalent to R^ values) accounted for by entering static, 
seasonal, and date-specific variables in the order listed in Table 1.3. A 
similar pattern of R^ through the crop year 1973“74 Is evident In both the 
Landsat and Tassel led Cap spectral measurement spaces. This temporal pattern 
for the static variables consists of a relatively low R^ on Date ^ 1 
(emergence), a high R^ in Date # 2 (jointing minus), a drop in R^ on date 
^ 3 Hointing plus), another R^ peak at date If 4 (turning minus’ , and drop 
in R^ on the last date. Seasonal variables accounted for littie variance 
(after the effects of long term or static variables were removed) in wheat 
spectral response over the year. Date-specific variables peaked in performance 
on date ff 3* 

Tables 1.7 and 1.8 list R^ values resulting when the seasonal variables 
are not included In the ordered (stepwise) regression procedure, in this 
case, most of the variance formerly attributable to seasonal variables is 
now accounted for by date-specific variables. The temporal pattern of R'^ 
values apoears similar to that described in Tables 1.5 and 1.6. 
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An explanation for the blmodal R temporal pattern obtained In this 
analysis must be developed by relating wheat canopy development to the 
signature predictor variables used In the regressions* Detailed canopy 
development Information over broad areas was not available; however, 
knowledge of general wheat development patterns for Kansas, known physiological 
relationships, as well as interpretation of full frame Landsat CIK transparencies, 
and review of variable-specific sensitivity analysis results (Table 1.9 ■* 

1.18) suggest the following hypothesis for the R* pattern. 

At emergence (date ^ !), the wheat "canopy" w{U be composed of Isolated 
shoots; consequently spectral s Ignature wl 1 1 be domlnatfid by soil surface 
reflectance. Soil reflectance will In turn be a function of surface horizon 
texture and moisture content. Reference to Tables 1.9 and I.IO confirms 
the Importance of soil reflectance generally ( dry soil albedo for bare soil) 
and texture (strongly related to soil available water holding capacity) In 
accounting for variance In spectral response. These same tables did not 
show the soli surface wetting measure (precipitation In the four days previous 
to Uandsat pass date - l.e. Four-day precipitation) as Important since 
no precipitation occurred In this period. 

The degree to which soil dominates reflectance will be a function of 
cultural practices and crop development. Favorable Fall growth environments 
will favor shoot development and reduce soil spectral influence. 1973 
planting season degree-days , one measure of growth environment, was shown 
to be one of the most Important factors (Table 1.9) accounting for variance 
In Landsat bands. Cultural practices Including planting date, seeding density, 
planting pattern, variety mix, and irr Igatlon/fertI 1 Izatlon application will 
also affect the percent 'Average of soil by wheat and degree of additional 
coverage at pass date due to canopy development. These cultural pr/>ctlce 
factors were poorly represented in this analysis. Only a very general measure, 
cultivated land percent , was Included. This cultivation practice variable 
was found s I gn i f leant . ~ At present, we ascribe the relatively low average 
total r 2 for Date ^ 1 to a lack of these cultural practice variables - - 
variables particularly Important In determining soil dominance of wheat 
signature In this biophase. 

Results In Tables 1.11 and 1.12 show long term (30 year average) growing 
season degree-days , cultivated l and percent , and soil water holding capacity 
account for most of the wheat "Spectral variance on the second date (jointing 
minus). Apparently, the degree of wheat canopy development could, at this 
blostage In this crop year, be expressed largely as a function of climate 
and soil type. In fact, an average of 87.5 percent of the variance in wheat 
field mean spectral response In Landsat bands and 7^-3 percent In the first 
two Tassel led Cap bands was accounted for by the static stratification 
variables as a whole. This situation represents most closely that originally 
hypothesized for signature extension stratification. 
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Table l,9i, WMS . (Crop Year 1973 - 7^1) 

Spectral Sensitivity Analysis (R^ Values) 


Dates October 20, 23, 29 
No. OF Segments: 9 
Total No. of Fields- 76 




L4 

L5 

1. 

Cultivated Pct. 

.124 

.160 

2. 

Crop Diversity 

.015 


3. 

Water Holding Cap, 

.076 

.091 

4. 

L.T. Grow, Season Deg. -Days 

.056 

.071 

5* 

1 .T. fiRDW. RPA«;nw pRpr.TP. 

art r « w f t a *m t a'^ w « a t t a aaa * t a 

no 

. \ jii . 



(.274) 

(.322) 

6. 

Grow. Season Deg, -Days 



7. 

Grow. Season Precip. 



8. 

January Temp, 



9. 

Plant. Season Deg, -Days 

.196 

.223 

10. 

Plant, Season Precip. 

(.196) 

(.223) 

11. 

Dry Soil Albedo 

.105 

’.095 

12. 

Four-Day Precip, 



13. 

Scan Angle 



14. 

Avg. Slope 



15. 

7/5 Ratio 

•.027 ‘ 

.036 



(.132) 

(.131) 

Tota 

,L 

.602 

.675 


17 

.058 

.112 

.119 

.196 

■ nno 

I uUo 

(.529) 


,0U 

(. 011 ) 

.027 


.177 

(.199) 

. 74 /^ 


OTAL 


Table l.lOi 


OP IS 
OF POOR QUALITY 

(Crop Year 1973 - Z'l) 

* 

Spectral Sensitivity Analysis (R^ Values) 


Date I October 2Qj 23 j 2A 
No, OF Segments.' 2 
Total No. of Fields; 29 



■ 

TCl 

TC2 

TC3 

TCi| 

1. 

Cultivated Pct. 

.175 

.5i|8 

,03i| 

.002 

2. 

Crop Diversity 





3. 

Water Holding Cap. 

• .122 

0 

,115 

.073 

It. 

L.T. Grow. Season Deg, -Days 

• 




5. 

L.T. Grow, Season Precip. 







(.297) 

(.5118) 

(.1A9) 

(.075) 

6. 

Grow. Season Deg. -Days 





7: 

Grow, Season Precip, 





8. 

January Temp. 





9. 

Plant, Season Deg. -Days 





in, 

Plant. Season Precip. 





11, 

Dry Soil Albedo 


• 



12, 

Four-Day Precip. 





13. 

Scan Angle 




• 

w. 

Avg, Slope 





15. 

7/5 Ratio 

0 • 

.390 

.372 

.251 



(.000) 

(.390) 

(.372) 

(.251) 

Total 

.297 

.938 

.522 

.325 
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Table Itlli 


KANS AS (Crop Year 1973 - 7A) 


"^ORIGINAU PAGE e 

Of POOR QUALVVY 


■ 9 

Spectral Sensitivity Analysis CR^ Values) 



Date; April 18, 20 
No, OF Segments; 

Total No. of Fields* 56 





L5 

L7 

1. 

Cultivated Pct, 

.287 

.250 

.153 

2. 

Crop Diversity 




3. 

Water Holding Cap. 

.221 

.226 

.241 

A. 

L.T. Grow, Season Deg, -Days 

.296 

.251 

.410 

5. 

L.T, Grow, Season Precip, 

.1^1 

.096 

.055 


»• 

( QIiP.'s 

f .R90^ 

f 

■% • W r 

6. 

Grow, Season Deg, -Days 

.009 

.010 


7. 

Grow, Season Precip, 




8. 

January Temp, 




9. 

Plant. Season Deg, -Days 




10. 

Plant. Season Precip. 






(.009) 

(.010) 


11. 

Drv Soil Albedo 



.003 

12. 

Four-Day Precip. 




13. 

Scan Angle 




W. 

Avg. Slope 




15. 

7/5. Ratio 

.•025 

.109 

.101 



(.025) 

(.109) 

(.104) 

Total 

.980 

.9^10 

.963 
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Table 1.12’; 


KANSAS (Crop Year 1973 - 7^1) 


(original page is 

OF POOR QUALITY 


Spectral Sensitivity Analysis (R^ Values) 


Date: April 18. 20 

No, OF Segments: 4 
Total No, of Fields : 55 




TCI 

TC2 

TC3 

TC9 

1. 

Cultivated Pct. 

.110 

.323 

.001 

.234 

2. 

Crop Diversity 





5. 

Water Holding Cap, 

• .316 

,121 

.344 

.226 

i|. 

L.T, Grow, Season Deg. -Days 

.337 

.080 

.846 

.373 

5. 

L.T. Grow. Season Precip. 

.015 

.184 

.036 

.129 



(.778) 

(.708) 

(.867) 

(.962) 

6, 

Grow, Season Deg. -Days 

.001 

.027 

• .002 

.006 

7; 

Grow. Season Precip, 





8. 

January Temp. 





9. 

Plant, Season Deg. -Days 





10. 

Plant. Season Precip. 







(.001) 

(.027). 

(.002) 

(.006) 

11. 

Dry Soil Albedo 

.005 

.005 

.005 


12. 

Four-Day Precip. 





13. 

Scan Angle 




• 

W. 

Avg, Slope 





15. 

7/5 Ratio 

. .029 

.216 

.048 

.019 



(.029) 

(.221) 

(.053) 

(.019) 

Total 

.808 

.955 

.923 

.987 


1.19 




iM 3 r)fW ?‘H 6 v;ever, by date § 3 (Jointing plus) the static variables no longer 
/Recount for‘a high percentage of the spectral variance (Tables 1.13 and 
l.l^i), The band 7 to band 5 ratio , a date-specific variable, Is dominant. 

Only crop d Ivers I ty , and to a lesser extent long term growing season 
peecipl tat I on , cultivated land percent , soil water holding capacity , and 
1 97^ g rowing season cleg ree^^^ account for any more than neg lib le amounts 

oT variance In wheat field mean response. At this time of year, the wheat 
canopy covered close to 100 percent of the soil surface In most western 
and central portions of Kansas. Canopy near Infrared reflectance (band 7 
In this case) was at Its peak or close to either side of the peak In most 
sample segments. A similar though reverse pattern (trough instead of 
peak) existed for chlorophyll -related depression of reflectance In band 5 . 

As a consequence, the average wheat spectral response curves 

appeared very flat when viewed over different climatic, land use, and soil 

association strata. This Is not to Imply that differences do not exist. 

On this date segments IO 36 and IO 65 (Southwest Crop Reporting District) 
have very similar Landsat band 7 to band 5 ratios, a measure used In this 
study as an Indicate' of crop development stage (maturity). Segment 1109 
Is, on the other hand, tracking spectrally behind 1036 and IO 65 with a 
higher 7 to 5 ratio. Thus differences are present, but not as pronounced 
as on dates 2 and § k. 

The spectral differences, or variance, on date # 3 that were explained 
by static stratification variables can be given the following Interpretation. 

As just mentioned, the largest spectral differences occurred between segment 
1036 and 1065 . Physically, 1109 differs from the former two by having 
denser texture soils with higher water holding capacities. This characteristic, 
when combined with higher growing season rainfall, leads to delay. In crop 
maturation relative to IO 36 and IO 65 - - a situation confirmed by the 
spectral data. In addition, the climatic environment In 1109 Is sufficiently 
different from 1036 and 1065 such that one would suspect that different 
varieties of wheat more appropriate to the humid environment would be planted. 
This may manifest Itself In wheat varieties with delayed maturation behavior. 

All three factors - growing season, precipitation, crop diversity, and 
available water holding capacity were found to account for variance in the 
analys Is on date § 3. 

The majority of the explained variance In the Landsat bands and In 
the Tassel led Cap band 2 was attributed to the Landsat band 7 to band 5 ratio. 
Importance of this date-specific variable relative to static stratification 
variables supports the conclusion that local growth environment variations 
rather than major differences In climate and soils dominate on third date. 

One further relationship should be noted for date § 3. At this time segment 
1109 was behind IO 36 and IO 65 In 197^ accumulated growing season degree-days, 
a situation favoring later wheat maturation In 1109. This seasonal degree- 
day variable was Identified as accounting for a small amount of variance, 
particularly in Landsat band 7 and Tassel led Cap band 1. 


if.- 

L3, KANSAS (Crop Year 1973 - 74) ofS’qwu,? 

^ Spectral Sensitivity Analysis (R^ Values) 

Date: *^ay 6, 9 


No. 

OF Segments; 3 




Total No, of Fields; 62 






L4 

L5 

17 

1. 

Cultivated Pct, 

.086 

.143 

.001 

2. 

Crop Diversity 

.151 

.129 

,023 

3. 

Water Holding Cap. 

.073 • 

.105 

.039 

4. 

L.T. Grow. Season Deg. -Days 




5. 

L.T. Grow, Season Precip. 

.004 

.011 

.110 



(.314) 

(.388) 

(.173) 

6. 

Grow. Season Deg. -Days 

.025 

.013 

.091 

7. 

Grow, Season Prec!?, 




8 . 

January Temp. 



• 

9. 

Plant, Season Deg. -Days 




10. 

Plant. Season Precip. 

> 





(.025) 

(.013) 

(.091) 

11. 

Dry Soil Albedo 

.001 


.002 

12. 

Four-Day Precip. 




13. 

Scan Angle 




14. 

Avg, Slope 




15. 

7/5 Ratio 

.342 

.364 

.257 



(.343) 

(.364) 

(.259) 

Total 

.681 

.766 

.522 
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Table 1,14! 



(Crop Year 1971’ - 74) 


Spectral Sensitivity Analysis (R^ Values) 


Date! May 6, 9 
No, OF Segments! 3 
Total No, of Fields! 62 




TCI 

TC2 

TC3 

1, 

Cultivated Pct. 

.030 

,101 

.453 

2, 

Crop Diversity 

.150 

,047 

.006 

3, 

Water Holding Cap, 

• .082 

,023 

,498 

4. 

L.T. Grow. Season Deg, "Days 




.5, 

l.T, Grow, Season Precip, 

.006 

,044 

,007 



(.267) 

(,215) 

(.964) 

6. 

Grow. Season Deg, -Days 

,09^ 


.005 

7, 

Grow, Season Precip, 




3, 

January Temp, 




9, 

Plant, Season Deg, -Days 




10, 

Plant, Season Precip. 






(,094) 


(.005) 

11, 

Dry Soil Albedo 


,002 


12, 

FouR"Day Precip. 




13, 

Scan Angle 




14, 

Avg, Slope 




15, 

7/5 Ratio 

,080 

,405 

.001 


1 

(,080) 

(.467) 

(.001) 

Total 

,442 

,682 

.670 


TC4 

,017 

,031 


(,048) 

,026 


(,26) 


,038 

(,038) 

,111 


msm 


r 




By dote ^ ^ (turning minus) the static variables have reasserted their 
Influence. Differences In soil water holding capacity combined with 
precipitation Input ( long term and seasonal precipitation) have apparently 
(see Tables 1J5 and l .Tb) dll^ferentlal ly stressed the wlVeat crop. This 
stress when combined with differing energy Input ( long term growing season 
degree-days) explained approximately one third to one hall* of the total 
spectral variance In all Landsat and Tassel led Cap bands. Most of the 
remaining explained variance on date # l| was associated with the landsat 
band 7 to bond 5 ratio. As such, this variable represented date-specific 
and location-specific crop maturity Influences (natural and management 
practice related) accounted for by static or seasonal variables. 

Usefulness of results from date # $ (turning plus) are limited by 
the fact that two segments with available data, IO 36 and IO 65 , fall In 
the same climatic stratum. Three variables account for most of the explained 
variance (see Table 1.17); In order of Importance these are crop diversity , 
Landsat band 7 to 5 ratio , and 197A growing season degree»days accumulated 
to the date In quest Ion ( Band I ng bn band 6 prevented genera t Ion of 
Tassel led Cop spectral values. 

Much of the explained difference can probaibly be attributed to the 
larger seasonal degree-day Input Into segment IO 36 relative to IO 65 , and 
also to date specific (indicated by 7/5 ratio) natural or cropping practice 
Influences at the field level, The Importance of crop diversity Indicated 
In Table 1,17 may be an artifact of the limited data set for this date. No 
physical explanation cat at present be given for the Importance of this 
variable. 


Table 1.18 summarizes the average (proportion of total spectral 
variance) values over all five dates for Landsat bands 5 and 7 and for 
the two. most Important Tassel led Cap bands I and 2. Landsat band A 
results were very similar to those for band 5 and were thus omitted. The 
grand average column on the right side of Table I.I 8 Indicates that 
cultivated area percent , soil available water holding capacity , long term 
growing season degree-days , and Landsat band 7 to "band' ~5 r^^ V'/ere tlie 
most Important signature pr ed 1 c tor va r I ab I es . Especl a 1 1 y notT worthy Is 
the fact that the band 7 to 5 ratio, In spite of being the last entered 
In the regression equations, accounted for the largest amount of spectral 
variances In all but the first Tassel led Cap band. 


Given the entry of variables In the order listed In Table 1.18, we 
can roughly express the results averaged over all dates as follows. 

Cultivated area percent, soil available water holding capacity, and long term 
growing season degree-days each account for 1/8 of the total spectral variance. 
Crop diversity and long term growing season precipitation each explain about 
1/16 of the variation. Landsat band 7 to 5 ratio accounts for approximately 
3/16 of the variance, while the remaining variables together explain another 
1/16. Altogether, an average of approximately 3/^ of the total spectral 
variance was explained by variables In the regression equations. 
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Table 1.15; 


KANSAS (Crop Year 1973 - 7A) 

o 

Spectral Sensitivity Analysis (R Values) 


Date ; May 2^1. 26. 27 
No, OF Segments; 6 
Total No. of Fields* 110 



• 

L4 

L5 

L7 

1. 

Cultivated Pct. 

,011 

.008 

.072 

2. 

Crop Diversity 

.012 

.042 

.078 

3. 

Water Holding Cap. 

.158 

.209 

.008 

1 A. 

L.T. Grow. Season Deg, -Days 

.312 

.296 

.086 

i 5. 

L.T. Grow. Season Precip, 

.018 

.032 

.280 

f 


f 

f t:QQ\ 
\ . 

f c:9/l^ 

\ 1 

1 6. 

Grow, Season Deg, -Days 

.010 

.011 

.008 

t 7. 

Grow. Season Precip. 

.053 

.064 

.001 

i 8. 

January Temp. 



, 

9. 

Plant, Season Deg, -Days 

.005 

.004 

.057 

10. 

Plant, Season Precip. 

.008 

.010 

.058 



(.076) 

(.089) 

(.124) 

11. 

Dry Soil Albedo 


• 


12. 

Four-Day Precip. 




13. 

Scan Angle • . 




14. 

Avg. Slope 




. 15. 

7/5 Ratio 

..168 

.152 

.088 

t 


(.168) 

(.152) 

(.088) 

Total 

.755 

.829 

.737 


jn 
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Table 1.16! KANSAS (Crop Year 1973 - 74) p°0R Quality 

Spectral Sensitivity Analysis (R^ Values) 


Date: May 24. 26. 27 

No. OF Segments: 6 
Total No. of Fields: 110 




TCI 

TC2 

TC3 

TC4 

1. 

Cultivated Pct. 

.032 

,018. 

.032 

.037 

2. 

Crop Diversity 

.004 

.130 

.008 

.066 

3. 

Water Holding Cap. 

■ .132 

.162 

.083 

.164 

^1. 

L.T, Grow. Season Deg, -Days 

.352 

.034 

.165 

.049 

5: 

LJi Grow, Season Precip. 

.002 

.298 

.057 

.223 



(.522) 

(.642) 

(.345) 

(.539) 

6. 

Grow, Season Deg. -Days 

.001 


.009 

.001 

7; 

Grow. Season Precip. 

.075 

.024 

.102 

.067 

8. 

January Temp. 





9. 

Plant. Season Deg. -Days 

.020 

.014 

.018 

.023 

10. 

Plant. Season Precip. . 

.033 

.005 

.048 




(.129) 

(.043). 

(.177) 

(.091) 

11. 

Dry Soil Albedo 


• 



12. 

Four-Day Precip. 

. 




13. 

Scan Angle 




• 

1^1. 

Avg, Slope 





15. 

7/5 Ratio 

.048 ■ 

.266 


.169 



(.048) . 

(.266) 


(.169) 


Total .700 .950 .522 .799 
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Table l.lZi 



(Crop Year 1973 


74) 


Spectral Sensitivity Analysis (R^ Values) 


Date: June 14 
No, OF Segments; 2 


Total No, of Fields* 47 






L4 

L5 

L7 

1, 

Cultivated Per. 

.001 


.010 

2. 

Crop Diversity 

.275 

.312 

.223 

3. 

Water Holding Cap. 

.012 

.010 

.064 

H , 

L.T, Grow, Season Deg. -Days 




5. 

L.T. Grow, Season Precip, 

( .‘288) 

/ rrnrw 

\ tOll) 

(.298) 

6. 

Grow, Season Deg. -Days 

.100 

.095 

.118 

7. 

Grow. Season Precip. 




8. 

January Temp. 




9. 

Plant. Season Deg, -Days 




10. 

Plant, Season Precip. 

(.100) 

(.095) 

(.118) 

11. 

Dry Soil Albedo 




12. 

Four-Day Precip. 




13. 

Scan Angle . . 




W. 

Avg, Slope 

• 



15. 

7/5 Ratio 

.111 

.158 

.003 


(.111) (.158) (.003) 

.500 .574 .419 


Total 
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Table 1.18; Summary of Ordered Contributions to R by Band Averaged 

Over Oates 
0 


Band ; Row 




L5 

L7 

TCI 

TC2 

Average 

1. 

Static Variables 







1 . CULTPCT 

.111 

.059 

.087 

.227 

.121 


2. COIVSITY 

.097 

.089 

.038 

.044 

.067 


3. AWC 

.128 

.099 

.163 

.074 

.116 


k. LTGSDD 

.124 

.139 

.172 

.029 

.116 


5. LTGSP 

.028 

.091 

.006 

.131 

.064 

11. 

1973“?^ Season- 
Specific Variables 







6. SUMGSDD 

.026 

.043 

.024 

.007 

.025 


7. SUMGSP 

.013 


.019 

.006 

,010 


8 . JANTEMP 

9. SUMPSDD 

,045 

.014 

.005 

.004 

.017 


10. SUMPSP 

.002 

.012 

,008 

.001 

.006 

111. 

Date-Specific 

Variables 







1 1 . ALBDRY 

,019 


.001 

.002 

.005 


12. PPT^DA 

13. SCANANG 

14. SLOPE* 

15. RASF 

.164 

.125 

.038 

.334 

.165 


* 


F 1 xed : 


Not date-specific but placed here for purpose of determining the 
2 

R overlap between scan angle and slope effects. 


m 






mmmmm 


WPP 






Examining the spectral predictor variable Importance by band, the 
-ol lowing relationships were evident. Landsat band'i 5 and 7 and Tassel led 
Cap bands 1 end 3 are again used as the most lmpv«*(«nt measures of spectral 
response. On the average, cultivated area percent , sol 1 aval lable water 
noldlng capacity , and long term growing season degree-days were the static 
»^arlables cons 1 s tan 1 1 y account I ng for sign i f 1 can t amounts of spectral variance 
in band 5. Much of the remaining explainved variability In band 5 was 
accounted for by the band 7 to band 5 ratio. Long term growing season 
degree-days dominated explained variance In Landsat band 7 on dates 1 and 2, 
while long term growing season precipitation Is the most significant static 
v-arlable on dates 3 and V. Landsat band 7 to band 5 ratio Is consistantly 
Important over dates 1 to k and dominates on dates 3 and 4. Date ^ 5 Is 
excluded from this analysis due to the low number of segments and poor 
distribution between climatic and soil types. 

Soil available water holding capacity and long term growing season 
degree-days are the most Important variables In explaining Tasselled Cap 
band I variance on dates 2 and 4. Cultivated area percent was also Important 
earlier In the growing season (dates 1 and 2). Crop diversity was significant 
on date 3 as well as on date 5. In the case of Tassel led. Cap band 2, the 
Landsat band 7 to band 5 ratio v;-as consistently an Important predictor 
variable. This result supports the hypothesis that the 7 to 5 ratio is 
strongly correlated to the Kauth gressn band. Even after the effects due to 
all other variables were removed, the Landsat 7 to band 5 ratio still accounted 
for 46.5 percent of the total variance In Tasselled Cap band 2 on date # 3. 

This time of year corresponded approximately to the spectral wheat calendar 
peak for the 7 to 5 ratio In western and central Kansas. Cultivated area 
percent was also Important In explaining variability in Tasselled Cap band 2 
on dates 1 through 3. On date ^ 4, crop diversity , soil water holding 
capacity , and long term growing season precipitation joined the Landsat band 
7 to band 5 ratio In accounting for spectral variance in the second Tasselled 
Cap band. 

Regression Without Prior Ordering 

A second regression analysis was performed to determine the most 
significant signature predictor variables when no_ order of entry was specified 
In advance. For this purpose, a forward selection procedure was applied to 
the data used In the previous ordered regression analysis. 

The product of this second analysis was a ranking of signature predictor 
variables from first entered (into the regression) to last along with the 
Increment in spectral variance explained (R^ Increment). Results were 
developed from separate regressions by band by date. A summary of the three 
most important variables at the top of each list is given in Table 1.19. 

The list together with corresponding R Increment and F-statlstIc data 
were then used to rank the signriture predictor variables In terms of relative 
Importance, highest to lowest. For the 15 predictor variables defined Iri 
Table 1.3 and also listed In Table 1.18, the following rankings resulted; 

1. Landsat band 7 to band 5 ratio (RASF) 

2. average bare soil albedo (AlBDRY) , accumulated precipitation 
for t~e four days previous to pass date (PPT4DA) , Landsat 
scan a^gle (SCANANG) , and slope percent (SLOPE) 


Table 1.19: Summary of Three Most Important Signature Predictor Variables by Band and by Date Using 

Forward Regression (No initial Order of Predictor Variable Entry Specified). Increment 
to Shown in Parenthesis. 
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3, long term growing season degree-days and precipitation (LTGSOD 
and LTGSP) 

cultivated area percent (CULTPCT) , crop diversity Index (CDIVSITY), 
and soli available water holding capacity (AWC) 

5. 197^ growing season degree-days and precipitation accumulated to 

the Landsat pass date (SUMGSDD and SUMGSP) , average January 197^ 
temperature (JANTEHP) , and Fall 1973 planting season degree-days 
and precipitation (SUMPSDO and SUMPSR) 

Brief comments regarding these rankings follow; 


1. BASF was clearly the most Important variable. It ranked on the 
top of the list for 11 of the 28 regressions under consideration; 
In 13 more regressions It ranked second. Roughly speaking, RASF 
explained approximately half of the spectral variance accounted 
for by predictor variables. Its dominance Is largely due to 
the fact that it was a scene and band dependent variable varying 
by pixel. All other variables, varied only by segment or 
stratum within segment. 


2. ALBDRY, PPT^DA, and SCANANG represent the remainder of the date- 
specific variables examined in this study. Together with SLOPE, 
these variables as a group exhibited erratic ranking behavior. 
However, all but SCANANG came out on top In at least one regression. 
They all tended to be very highly correlated with each other and 
with several members of groups 3 and 5. SCANANG seemed to be the 
least important. Their appearance., even at the top of a list, 
seemed highly dependent on the presence or absence in the list 
of other variables with which they were highly correlated. 


3. The static variables LTGSDD and LTGSP exhibited behavior similar 
to group 2, except they gave a high negative correlation (r = -.9) 
for date 4, which contained 6 of the 7 segments. They each came 
out on top in at least one regression. LTGSDD appeared slightly 
more important, and together they were slightly less important 
than group 2. 


4. CULTPCT, CDIVSITY, and mWC were consistently of medium importance. 
These variables tended to be highly correlated with each other 
and only moderately correlated with other variables. 


5. The seasonal variables as a group generally received lower rankings. 
Notable exceptions to this trend occurred in band TC 3 on dates 
4 and 5. Interestingly, these variables did somewhat better in 
Tassel led Cap as opposed to Landsat spectral space. 
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SECTION 2.0: 

DEVELOPMENT OF PHOTO I NTERPRETATI ON METHODS BASED 


ON MULTI TEMPORAL LANDS AT DATA 


2.0 


DEVELOPMENT OF PHOTO INTERPRETATION METHODS BASED ON MULTITEMPORAL 
LAND SAT DATA 


2.1 summary task description 

« 

The objective of Task II Is to develop photolnterpret Ive rethods 
for utilizing mu 1 tl temporal Landsat data which are Improvements to the 
current Lade methods. In order to accomplish this objective, Task II 
has been approached through three subtasks. These subtasks are: (1) Subtask 
A - Familiarization with JSC/LACIE Photointerpretation Procedures, 

(2) Subtask B - the Development of Mul tl temporal Interpretation Procedures 
whereby Individual temporal Images and temporal spectral data are analysed 
by means of a decision logic for the Identification of wheats and (3) Subtask 
C - the Development and Evaluation of Methods for Reducing Multitemporal 
Data to a single Image set of mul tl tempora 1 ly combined/enhanced Images and/or 
a numeric representation of the spectral data. 

2.2 SUBTASK A; FAMILIARIZATION WITH JSC/LACIE PHOTO INTERPRETATION PROCEDURES 
Objective 

The objective of this subtask is to familiarize UCB with JSC/LACIE 
photointerpretation procedures. In addition to an initial familiarization 
effort with current LAC IE procedures, a continuing effort Is maintained to 
remain current on all Implemented and proposed photointerpretation procedure 
modifications. This subtask Is vital to efficient and effective performance 
upon the other two subtasks within Task II. 

Approach 

An Initial familiarization visit to JSC was made In August 1976. 

During this visit CAMS Operations personnel provided tutorial sessions and 
over-the-shoulder-look interpretation sessions for the RSRP personnel. In 
addition to reviewing CAMS procedure documents as they become available, 

UCB will maintain contact with CAMS Operations personnel and will request 
additional Interactive sessions as needed to remain current on any modifications 
to the JSC/LACIE photointerpretation procedures. 

Progress Last Quarter 

A member of the RSRP staff (C.M. Hay) spent three weeks (28 March 1977 
to 15 April 1977) at JSC aiding in the development of the LIST method. This 
provided an excellent opportunity for contact with LAC IE operational analysts 
and a better appreciation of the ana1yst‘s interpretation environment problems 
was acquired. Also an update of progress of new Interpretation procedures 
(Small Fields Procedure and Procedure 1) was gained. Contacts such as these 
with the operational environment of LACIE are the most efficient method for 
UCB personnel to stay current with JSC/LACIE photointerpretation procedures. 


m 




During this quarter the special PFC Alternative Product Evaluation 
was concluded. This secondary study provided the RSRP an opportunity to 
compare the usefulness of current and potential future Interpretation 
products to the analysts. As a result of this study an increased awareness 
of the significance of data distortion as occurs when Landsat spectral 
data Is mapped to display or Image production was gained and subsequently 
will benefit the other two subtasks of this task. The results of the PFC 
Alternative Product Evaluation are presented In Appendix A. 

Work To Be Accomplished Next Quarter 

Procedure 1 Is scheduled for Implementation within the next quarter at 
JSC. After implementation has occurred, UCB plans a visit to JSC for a 
briefing from operational analysts and other JSC personnel on Procedure I 
as actually Implemented. Problems pertaining to analyst labeling that may 
be associated with Procedure I will be reviewed at this time. 

2,3 SUBTASK B; DEVELOPMENT OF MULTITEMPORAL INTERPRETATION PROCEDURES, AND 

SUBTASK C: DEVELOPMENT OF IMAGE REPRESENTATION TECHNIQUES FOR COMPRESSED 

MULT I TEMPORAL DATA. 


Background 


In order to develop products and procedures that can potentially Increase 
the labeling accuracy of analysts working with Landsat data such as in LACIE, 

It Is first necessary to understand the components of the labeling process. 

The labeling process can be considered to consist of two main components 
^sature detection and (2) feature Identification . Feature det(»,ctlon 
Is deft ned to be the action of dTscriml nation, based on spectral, spaclal 
and temporal characteristics, a unique entity observable within the Landsat 
multitemporal spectral data. Feature Identification Is defined to be the 
action of assigning a name (eg. wheat, non-wheat) to the' detected feature, 

Correct feature Identification can not properly proceed unless feature 
detection has first occurred. Feature detection, however, does not Insure 
feature Identification. Errors In labeling can thus occur due to (1) failure 
to detect a feature of Interest, and (2) failure to correctly Identify a 
detected feature. 

Errors Due to Failure to Detect a Feature of Interest 

In multitemporal analysis of Landsat (low resolution) data for crop 
Identification, the characteristics that the analyst seeks to detect are (1) the 
presence or absence of a vegetation canopy within specific biophases, (2) the 
relative condition or stage of development (quality of presence) of a 
vegetation canopy within specific blophases, and (3) the spaclal distribution 
(local field size and shape, and overall proportion of similar fields In the 
segment) of a given spectral/vegetation canopy type. Of these three characteristics 
the first Is the most Important for wheat detection and Identification. 
Determination of the other two characteristics Is necessary when significant 
overlap between wheat and confusion conditions exists or when acquisitions 
are missing or of poor quality. If a vegetation canopy can not be detected 


w!thln a critical wheat blophase the probability of ultimately Identifying 
the feature as wheat Is lowered. Thus If the standard PFC Product 1 fails 
to represent a vegetation canopy In a manner normally expected by the analyst, 
there Is a danger of mislabeling due to failure to detect the presence of 
vegetation on a given acquisition. Product 1 tends to Inadequately represent 
low density vegetation canopies. Many> though not all, of these low canopy 
situations are detectable upon examination of the actual numeric spectral 
data. Thus Auxiliary products which are designed to more clearly Indicate 
the presence of low density canopies or other difficult canopy states should 
aid In Increasing labeling accuracies where errors are due to non-detection. 

In assessing the quality of the vegetation canopy, the numeric spectral 
data allows the analyst to differentiate more finely between stages of canopy 
development or condition than does the standard Product 1. This capability 
may be useful In discriminating wheat from close confusion crops. Auxiliary 
products that allow the analyst to evaluate the actual numeric spectral data 
Improve detection and measurement of canopy presence or absence as well as 
canopy condition. 


Products currently being Investigated and refined as auxiliary aids to 
the analyst are listed in Table 2.1. The numeric spectral data products are 
to give tne eneiydi. a picuisic; mcasuic? wi s.»v./fiv4 1 u i wn anw umi 

callbra’te the standard Images with which he Is working. The Image display 
products are meant to give the analyst a convenient representation of the 
spaclal distribution of a given condition. 


Errors Due to Incorrect Identification of a Detected' Feature 

While the standard Image and auxiliary products allow the analyst to 
detect a feature, It Is actually ancillary data from outside a specific 
set of Landsat spectral data that allows the analyst to Identify a detected 
feature. The analyst must make correlations between a given spectral -temporal 
response pattern and the ancillary data about ground conditions. An analyst's 
ability to do this Is dependent upon his educational background, his experience 
with the environment that he must analyze, the quality of training, the 
definitiveness of data correlations that have been given to him as guidelines, 
and the quality of his own analytical thought processes. At present, the 
analyst has very little data with which to make the needed spectral to 
ground correlations. Landsat data has not been available for a period 
sufficient to develop all the necessary correlations. The analyst, therefore, 
must develop these correlations (or appropriately inferences to possible 

correlations) for himself from the ancillary data and spectral data given to 
him. Until the necessary spectral to ground condition correlations can be 
determined and def ini tl'vely presented to the analyst the Identification 
process will continue to be heavily dependent on the skills and experience 
of the Individual analyst. 
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Table 2.1; Data To Be Evaluated for the DEVELOPME^^• of Products for 
Improved Feature Detection 

’ Numeric Spectral Data 

• Spectpal Means^ Standard Deviation and Covariance 

• 7/5 Field/Pixel-Specific Ratios 

• Kauth Green Numbers and Brightness Numbers 
• Image Display 

' Cluster Maps 

• Principal Components Representation 
' Kauth Rotated Data 

• Texas A & M University Data Reduction Technique 
■ 7/5 Ratio 

• Multitemporal Combinations of Above 
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The ancillary data currently considered to be the most valuable to the 
analyst In aiding him to develop the needed correlations are listed In Table 
2,2. This ancillary data should not only contain mean values and descriptions 
of average normal conditions but also should Include variability Information 
about the data. The variability data should Include temporal variability 
(year to year) as well as spaclal variability if at all possible. 

2.3.1 SUBTASK B: DEVELOPMENT OF MULTITEMPORAL INTERPRETATION PROCEDURES 

Objective 

The purpose of this subtask Is to develop Improved mu 1 1 1 temporal Interpre- 
tation procedures which will increase Al labeling accuracies relative to current 
LAC IE Interpretation procedures. Specifically procedures which allow the 
analyst ‘to utilize both standard mul tl tempora I image presentations and 
preprocessed spectral data presentations are being explored. 

Approach 

This subtask builds upon Landsat Image Interpretation procedures 
developed In a previous UCB-LACIE task (see Hay and Thomas, 1976) and 
endeavors to Incorporate the use of more detailed spectral Information 
(such as products being developed in Subtask C) along with the analysis 
of the standard Image products. The time has come for more frequent 
usage of the detailed spectral information actually contained In Landsat data. 
Standard Image products such as Product I do not allow for the refined 
feature discrimination that Is possible with Landsat data. This is not to 
say that standard image products such as Product 1 are unnecessary; they 
are vital. What Is now needed, however. Is the Incorporation o^^ more 
refined spectral data into the analysis procedures along with standard 
Image products. 

A fair amount of additional Information Is available In the actual 
Landsat spectral data. However, formats for effective presentation of 
this data and methods of Incorporating It Into the analysis procedure need 
to be developed. Subtask C Is addressing the development of effective 
presentation formats. 

Subtask B addresses the development of methods and guidelines for 
the Integration of more refined spectral data Into the analysis procedure. 

A problem In developing these guidelines is the lack of sufficient 
correlation of Landsat spectral responses to given ground conditions. 

Spectral response-to-ground condition relationships are just now beginning 
to be established. As yet, the results are meager and inconclusive. This 
should not stop the use of more refined spectral data by the analysts 
however. Some general guidelines can be developed now for the analysts 
and refined later as more definitive results from data correlation studies 
become aval lable. 


Table 2,2 i Ancillary Data flosT Necessary for Feature Identification 

, Data for all Crop IIade Available to AI 

• Crop Calendars 

' Crop Historical Proportions 
' Cropping Practices 

' Meteorological Data Affecting Crop Development and/or Spectral Response 
OF Crop 

• Recent Precipitation 

• DnxcrMTTAi Yte! ri 
I VI uni mu i icLu 

' Episodal Events Data 

Spectral Response Correlated to Crop Development Stages 



Procedures are being explored that make use of vegetation Indication 
measures such as the band 7 to band 5 ratio or Tassel led Cap green and 
brightness numbers, information about the temporal pattern of the 7/5 ratio, 
green numbers, and brightness numbers for wheat and other crops as well as 
thresholds for the presence of vegetation can be utilized by the analysts. 
This task seeks to develop procedures and guidelines for the effective 
Incorporation of this data Into the analysis procedure. 


Another aspect of the labeling procedure being explored In this 
subtask Is that of sampling within the segment. The development of 
procedures for better data usage and more effective sampling are being 
pursued In full recognition of the development of Procedure 1 and with 
the view that any procedures developed would, theoretically, be compatible 
with the philosophy behind Procedure 1, 


All procedure modifications will be evaluated against the JSC/LACIE 
field Interpretation procedures (the control treatment). Evaluation criteria 
for the procedure modifications Include(l) ability to Improve labeling 
performance, (2) ability to Improve segment wheat estimate, (3) ability to 
decrease effort expended and Increase throughput rate, and (4) the 
repeatability of the results, These evaluation criteria are summarized In 
Table, 2.3. 


Table 2.3 

Subtask B: Development of Mu 1 1 1 tempera I Interpretation Procedures 


Evaluation Criteria 

Analyst Identification Accuracy 
Wheat Estimate for the Segment 
Effort Expended /Throughput Rate 
Repeatability of Results 


Progress Last Quarter 

Ground data for the 1975-1976 Kansas and North Dakota segments was 
received, catalogued, and filed. Ten (10) segments (5 from Kansas and 5 
from North Dakota) from the 1976 data set that will be used for procedure 
evaluation were processed this quarter by the JSC/LAC! E fields procedure. 


2.7 



1. Define spectral subclasses by analysis of the full segment. 

2. Select training fields from within spectrally homogeneous 
areas to represent all subclasses defined in Step 1. 

3. Identify all spectral subclasses as wheat or non-wheat 
jslng multi-temporal Interpretation procedures. 

4. Select five test fields which have not been selected as 
training areas. 

5. Digitize and verify training and test field coordinates 
using a coord Inatograph and the LARS terminal. 

6. Submit segment for batch classification processing. 

7. Evaluate the classification results using class map and 
performance matrix for training and test areas. 

8. If necessary, modify training and submit for reclassification. 

To simulate the JSC/LACIE control procedures, UCB analysts employed 
the following procedure. 

UCB - simulation of JSC/LACIE procedures: 

1. On the RFC transparencies, define subclasses according to 
JSC procedures. 

2. Select training fields according to JSC procedures. 

3. Identify as wheat or non-wheat all subclasses using JSC 
procedures . 

4. Systematically select fifty test areas throughout the 
segment. (It is UCB's view that the five test areas 
as required by JSC procedures are insufficient for an 
adequate evaluation of the classification results.) 

5. Extract training and test field coordinates using 
the UCB-Remedys color monitor display system. A 
coord inatograph as used at JSC Is not available at UCB. 

6. Submit the training deck to the supervised classifier 
(CALSCAN) for processing. 

7. Evaluate the classification results using the class map 
and performance matrix for training and test areas. 

Results for the fields procedure, which will be used as the control treatment 
for later comparison with alternative procedures, are presented in Table 2.4. 

These results Include analyst labeling accuracies and proportion estimates, as 
well as accuracies and proportion estimates for CALSCAN* runs using Ai labeled 
fields as training (JSC fields procedure). 



2.8 


ABLE 2.^; Interpretation Accuracies Using JSC/LACIE Phase II Fields Procedures (Control) 
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Work wa^ begun t!hi 5 ‘quarter on the development of a procedure whereby 
the analyst uses output from mu I tl temporal clustering to produce a stratification 
of the data Into wheat probability strata. This procedure, called the Del ta 
Function Stratification (DFS) Procedure , requires the A I to analyze Landsat 
band 7 to band 5 cluster mean ratios from each acquisition date Included In 
each mul 1 1 temporal cluster mean vector. Based on the analysis the analyst then 
(1) groups the clusters according to the probability that a given mul t 1 tempora I 
cluster Is wheat and (2) defines a display function for the clustered data 
which will allow visual grouping of the above grouped clusters. Visual 
grouping of similar clusters within the display product facilitates the usefulness 
of clustered data to the analyst in the Interpretation process. The Al can more 
quickly analyze the data and relationships within the data will be more apparent. 

By stratifying the segment Into wheat presence probability strata, a 
separate estimate of mean and variance can be produced for each stratum thus 
potentially allowing more precise estimates for wheat to be made at the segment 
1 evel . 


Initially the Delta Function Stratification Procedure is being developed 
using a full season's set of acquisitions (i.e. at harvest procedure). However, 
once developed for an at harvest mode, continuing work will deal with application 
of the procedure to the production of mid and/or early season estimates. 

Basically the procedure Is as follows: 

1. ) The analyst selects for a given segment a set of 2 to 4 

acquisitions which he has determined to be both necessary and 
sufficient for the Identification of wheat. 

2. ) The mu 1 1 1 tempora 1 data is then clustered using the following 

procedure. 

a. ) Five 20 point by 20 line seed areas 

distributed throughout the segment* 
are clustered** for 10 iterations 
using a STDMAX » A.O, MAXCLS = 50, 
and NMIN » 30. 

b. ) Punched output statistics from the 

seed run are input into a second 
clustering run on the full segment for ' 
an additional 10 iterations at a 
STDMAX “4.0, MAXCLS = 60, and 
NMIN “ 30* This equates to a total 
of 20 iterations for the segment. 

c. ) Punched output statistics and a 

display map on tape are acquired 
from the second clustering run. 


The upper left corner of the five 20 point x 2 line seed areas are: 
a.) 1,1; b.) 1,97; c.) S8,48; d.) 176,1; 3.) 176,97. 


Clustering algorithm ISOCLAS - UCB's adaptation of JSC ISOCLS, 
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The stat deck from the second clustering run Is Input to a 
program called CLODER* CLODER computes 2 times the ratio 
of the bad 7 to band 5 means for each cluster for each 
acquisition date submitted to the mult I temporal clustering 
algorithm. 

The analyst then selects a cluster ordering reference 
acquisition date . This acquisition Is the one on which 
the analyst determines that wheat Is most discrimlnable 
from most other cond 1 t lons/crop types present In the 
segment. 

Clusters are then ordered according to the 2 x band 7 to 
band 5 ratio of the cluster means. The highest ratio Is 
listed first, and the lowest listed last. The higher ratios 
correspond to the clusters which represent actively 
metabolizing green vegetation. 

After the mul tl temporal clusters have been ordered by the 
7/5 ratios on the reference date, the ratio of 2 x band 
7 to band 5 means for all other acquisitions that were 
processed are recorded In proper temporal order along 
side the reference date ratio for each cluster. 

The difference or delta between the ratios for each pair 
of adjacent acquisitions Is then computed (See Table 2.5). 

By knowing the temporal (biostage) function of the 7/5 ratio 
for wheat (see Figure 2.1) and all other major crops grown 
within a region, the analyst can then analyze the deltas 
between each processed acquisition and assign each cluster 
to a probability of wheat stratum. Three main strata are 
currently being used in the procedure. An option is available 
to the analyst to determine 2 substrata for each stratum 
making six the maximum possible number of substrata. 

The three main strata are: 

a. ) High probability wheat/small grains stratum: 

probability ^ 50 ^ that the cluster Is wheat/ 
small grains. 

b. ) Medium probability wheat/small grains stratum: 

^ 25 ^ to < 50 ^ probability that the cluster 
is wheat/small grains. This stratum seems to 
contain primarily pasture, alfalfa, and range 
clusters in winter wheat areas. 

c. ) Low probabi 1 i ty wheat/smal 1 grains stratum: 

< 25 ^ probability that a cluster Is wheat/ 
small grains. This stratum In a winter wheat 
area usually contains summer crops such as 
corn or sorghum or fallow conditions. 

Once the clusters are assigned to strata the analyst then defines 
the display function for the clusters. This function is defined 
.so that clusters within the same stratum or substratum will 
visually group together when displayed. This is necessary so 
that the analyst can efficiently analyze the spacial distribution 
of related clusters. Visual grouping can be accomplished in 
two ways: (1) by assigning the same color code to all clusters 
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Figure 2.1: Example of Band 7 to Band 5 Ratio of Winter Wheat Cluster Means through a Growing Season 

Compared to the Crop Calendar Plot 






*>"« , , 

In the same substratum (this does not optimally allow for 
analysis of Individual clusters ) or (2) by assigning adjacent 
colors (e.g. red and orange, yellow and yellow green, blue green 
and blue) to adjacent clusters within a substratum as ordered 
by their 7/5 ratios on the reference date. Figure 2.2 Illustrates 
display procedure 1 as described above and Figure 2.3 Illustrates 
display procedure 2. In Figure 2.2 all the clusters belonging to 
one substratum have been Individually displayed, clusters from 
the other substrata being displayed as a single color for that 
substratum only. This Is due to a current limitation in the 
RSRP IGOR device whereby approximately 25 unique, visually 
separable colors are available. It is planned to remove this 
limitation during next quarter's effort by expanding the number 
of visually separable colors to between 60 and 100. At present. 

It Is possible to work around the available color limitation by 
using multiple Images and displaying individual clusters for a 
given substratum on each image. 

Work To Be Pursued Next Quarter 

The Delta Function Stratification Procedure will continue to be developed 
next quarter. An expanded number of segments will be processed and analyzed 
using this procedure to determine If it will be workable for the full variety 
of situations represented by LACIE sample segments. In addition, refinements 
to the quantitative test design will be considered as the Delta Function 
Stratification Procedure continues to develop. 

2.3.2 SUBTASK C: DEVELOPMENT OF IMAGE REPRESENTATION TECHNIQUES FOR COMPRESSED 

MULTI TEMPORAL DATA 

Objective 

The purpose of this subtask Is to develop several alternative Landsat 
data reduction/compression spectral representations and to evaluate their 
role in mul tl temporal training procedures. These spectral representations 
function as additional sources of information to the analyst interpreter (Al). 

As such they should enable the analyst to better discriminate crop spectral 
density function components. If the spectrally discrlmlnable classes can then 
be correlated to ancillary data (crop calendar data, historical agricultural 
statistical data, etc.) improved labeling should be possible. 


Approach 

Several mul ti temporal spectral representations developed from Landsat 
digital data will be presented to the image analyst. Three of these will recel 
major emphasis in quantitative testing. They include*; (1) Landsat band 7 
to band 5 ratio multidate sequences, (2) Tassel led Cap (Kauth) green band 
multidate sequences, and (3) principal component representations of Tassel led 
Cap multidate sequences. 


*A11 Landsat data will be first corrected to reference sun angle and haze 
condition using ERIM XSTAR procedure. 
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Figure 2.2: Display Procedure 1 - Delta Function Stratification and Display 

Procedure was utilized to group the ^7 mu 1 1 1 tempora 1 clusters 
Into five wheat probability strata. 

High probability wheat straturp (H) , represented by red-violet areas, includes 
clusters 2 x 7 / 5 j^gy 5 ^ 1*5 and 2 x 7 / 5 j^j|y g ~ 

Medium probability wheat stratum I (Ma) , represented by red areas, Includes 
clusters 2x7/5j^gy g ^2.5 and 2x7/9jy|y 3 > 1.1; most probably alfalfa. 

Medium probability wheat stratum 2 (Mb), represented by blue areas, Includes 

clusters 2.5 > 2x7/5f^gy 6 ^ ^^^^^Julv 8 ^ probably alfalfa, 

pasture, range, or urban? ^ 

Low probability wheat stratum 1 (La), represented by green areas, includes 
clusters 2x7/5j^^g g £ 1*1 ^^d 2x7/5j|,iy 3 > l.l» most probably corn and 

sorghum. ^ ^ 

Low probability wheat stratum 2 (Lb), represented by white and black areas, includes 
c.lusters 2x7/5j^gy 6 :i !•! ^^^'^^Julv 8 ^ probably fallow and 

water. ^ y -* 


ground data 
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Sample of 209 grid points 


Figure 2.3: Display Procedure 2 - The Delta Function Stratification and 

Display Procedure was utilized to group 57 multi temporal clusters 
Into six small grains probability strata. 

Ha - high probability small grains substratum 1 
Hb - high probability small grains substratum 2 
Ma ® medium probability small grains substratum 1 
Mb - medium probability small grains substratum 2 
La - low probability small grains substratum 1 
Lb - low probability small grains substratum 2 

Strata Ma, Mb, La, and Lb have been grouped Into one color (brown) for display 
so that an adequate number of colors would be available for display of Individual 
clusters within strata Ha and Hb. Warm colors are assigned to the clusters that 
have the higher 7/5 ratios progressively toward the cooler colors with decreasing 
7/5 ratios. This allows visual grouping of similar clusters (similar 7/5 ratios) 
so that discrimination of landscape features Is preserved thus aiding the analyst 
In his analysis of the data. 
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Bet-'i'een date differences in Landscit 7/5 ratios or in Tassel led Cap green 
numbers will also be included as appropriate In (1) and (2) respectively, 

Formats of image representation to the analyst will include hardcopy 
prints, color monitor display, and/or numerical representations. Color 
hardcopy Is generated by the RSRP light emitting diode film annotator 
(160R). Television monitor display consists of single Images (one or more 
bands) or simultaneous display of up to four separate images (each with one 
or r.ore bands). 

An analysis of variance design will be used to quantitatively evaluate 
the effects of the three major ipectral representation types (defined above) 
on At labeling accuracy. To implement this design, a series of interpretation 
tests will be performed on 1976 LACIE blind sites in Kansas and North Dakota 
for which ground data is available. Analysts will be asked to first interpret 
the 209 grid intersections per segment using only the standard PFC film 
product and standard ancillary data. Next, segment grid intersections will 
be interpreted without standard products using only one of the three alternative 
spectral representations at a time. Each analyst will be randomly assigned 
one alternative product (multidate 7/5, multidate green number, or principal 
component of multidate green number) per segment to use for interpretation. 
Finally, alternative products will be combined with standard products and 
used to label the grid intersections for each segment. Analysts will interpret 
each segment once in this manner using the data combination containing the 
alternative spectral representation assigned to them previously. 

The test performance measure is defined as the change in labeling 
accuracy obtained with the al terciative- spectral representations or with the 
combination of alternative and standard products versus use of the standard 
products alone. These labeling accuracy changes by segment by treatment are 
the observations entered Into the analysis of variance. Thus the experiment 
is seen to consist of three treatment groups ; (1) standard product versus 
standard product (the control treatment; expected differece equals zero), 

(2) alternative product versus standard product, and (3) combination of 
alternative and standard products versus standard product. Relative variation 
in accuracy differences between treatments (reported by segment) is expected 
to be less than corresponding variation in individual accuracy values. 
Consequently, the treatment-difference approach is expected to give the most 
powerful test for detecting labeling accuracy improvements due to alternative 
spectral aids at a given level of statistical significance. 

Other spectral represenat ions will be evaluated qualitatively relative 
to potential labeling accuracy improvement by experienced image analysts. 

These additional spectral products include the Texas A & M mu 1 1 i tempora 1 
compression image as well as variations of the Landsat 7/5 ratio and Tassel led 
Cap space data. 
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Progress Last Quarter 

Work this past quarter focused on (1) completion of major aspects 
of the experimental design described above, (2) Implementation of the EfUM 
XSTAR sun angle and haze correction algorithms, and (3) development of 
Interactive color monitor display software to rapidly create and manipulate 
spectral enhancements. In addition, two LAC IE blind site segments were 
selected in Kansas (1035 and 104l) and two In North Dakota (1618 and 1633) 
for use in defining the specific data combinations and display formats for 
each of the three major spectral representations to be analyzed In the 
analysis of variance layout. 

Next Quarter 

A total of 17 test segments In Kansas and North Dakota will be 
interpreted by analysts using the standard PFC and ancillary data products. 
The resulting labeling accuracies will be used as the reference- performance 
measure In the quantitative evaluation of the Landsat 7/5 ratio and green 
number spectral representations. Formats for presentation of spectral aid 
materials used will be finalized. All test segments will be preprocessed 
by XSTAR procedures before multi date spectral enhancements are generated. 
Interpretation using the resulting spectral enhancements will be initiated. 
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By; Andrew S, Benson 
Intnvductfon 

Presently, the training data Input to the CAMS are wheat and nonwheat 
subclasses which the Al's Identify on a standard image type (Product 1) 
produced by the production film converter (PFC) . There is evidence, however, 
that the parameters Input to the PFC, viz., bias, gain, and spectral limits, 
may oe producing a biased representation of the spectral data on the resulting 
Imagery. It Is conceivable that (1) indications of early canopy cover present 
In the spectral data would not be manifested in the Imagery, (2) some of the 
wheat and nonwheat subclasses selected by the Al could have no correlation 
with spectral classes as seen by the classifier, and (3) adjoining segments 
from one Landsat frame would be produced using different parameters, thus 
producing varying representations of similar ground conditions. 

In an attempt to Improve the standard PFC image, three additional Image 
types (hereafter referred to as Kraus, Hocutt, Kaneko) have been produced on 
the PFC using different formulae to compute the input parameters of gain and 
bias. These formulas plus tha one used to compute Product 1 are given in 
Table A. 1. Of the four Image types, Product 1 is the only one in which the 

spectral ranges for each band are computed Independently of the others. With 
the Kaneko product, the spectral ranges for each band are computed separately 
based on the dispersion from the grand spectral mean for all bands. With 
Hocutt, one spectral range Is computed for all bands based on the maximum 
dispersion exhibited by any of the three bands. With Kraus, the complete 
range of Landsat data are used to produce the Image product. 

In January 1977, personnel from RSRP began an additional secondary task 
at the request of JSC technical monitors which was designed to compare the 
usefulness of the current PFC Image product with these alternative Image 
products. The following paragraphs document the results of this study. 

Objective 

The primary objective of this study was to evaluate the relative usefulness 
of Product 1, Kaneko, Hocutt, and Kraus image products for detecting and 
Identifying the presence of wheat. A secondary objective was to determine 
which of the three alternative products would be the best supplementary aid 
to Product 1 for the detection and Identification of wheat during the early 
growing season when wheat canopy density is low. 


Table 11 i^ormula^ Used to Compute Input Parameters to PFC Image Types* 


Product 1 


Kaneko 


Hocutt 


Kraus 


PFC A.* 


256 


Gain 


1 U.-L, 


A. = 


256 


I U.-L, 




256 




A = 

0 Pq 


PFC B,= (-1)L,A, 
Bias ' ' ' 


B| = (-l)LjA| Bq = (-I)LqAq None 


U,= P, + 3S, 


^1 " Pq ^1 *^o ^P| 


L, = Pi - 3S, 


L| - Pq - M, Lq = Min (P, - 3S,) 


M, = Max [Pj + 3S, - Pq, 


Pq - P| 3S,1 


. ^ *^2 ^ P3 


* Internal document from LACIE Project Review Session, January 25-28, 1977, 
NASA Johnson Space Center, Houston, Texas. 


Initially, the four PFC Image products were ranked qualitatively for' 
each date throughout the growing season using the following seven photo- 
graphic image quality criteria: 

1. Hue range (widest * 1, lowest “ ^) , 

2. Overall image contrast range (highest « I, lowest *» 4), 

3. Overall image density (lowest = 1, highest « A), 

k. Wheat signature intensity (highest *■ 1, lowest = 4), 

5. Hue saturation (highest = 1, lowest * 4), 

6. Ability to define field location — feature specific contrast 
• (best * 1 , poorest = 4) , and 

7. Ability to represent spectral classes — spectral data 
representation (best “ 1, poorest = 4). 

The first three criteria were applied to overall segment appearance, 
and the last four were applied to specific wheat and nonwheat fields for 
class evaluation within the segment. 

A quantitative evaluation of all acquisition dates from each segment 
was as follows: 

l. Histograms of the spectral data from bands 4, 5, and 7 for 
the entire segment were produced, 

2. Statistics were computed for each band, viz,, mean, standard 
deviation, 'and covariance matrix, 

3. The spectral limits which applied to each image type as 
outlined In Table A.l were computed, 

4. From a sample of the wheat fields given in the ground data 
the average number of wheat pixels falling outside the 
spectral limits for each image type were determined, and 

5. The spectral values of "anomalous appearing" pixels within 
wheat and nonwheat fields which were identified during the 
qualitative evaluation phase of the study were quantified. 

Approach 

In order to evaluate the relative merits of the standard PFC image 
products vs the three alternative image products, ground data, the four 
image products per date, and spectral data were examined for six Great 
Plains winter wheat segments and four spring wheat segments. A list of 
the segments and the acquisition dates available to the RSRP are given 
In Table A. 2. 
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Results 


Qualitative Evaluation 

The ranking of nine segments based on the photographic parameters 
Is given in Table A. 3. (Segment 1622 was excluded from this evaluation 
because no ground data were available). The performance index is based 
on. the average ranking of each product by segment computed from all 
acquisition dates. The ranking index is based on the simple numerical 
ranking of the performance Index for each segment. 

As is clear from Table A. 3, Product 1 was the best photographic product 
followed by Hocutt, Kaneko, and Kraus, respectively. This was to be expected 
because the procedure by which the gains and bias are calculated for Product 1 
will maximize the contrast and brightness of an image. In comparison, for 
many of the acquisition dates, the alternative products were poorly exposed, 
had poor color balance, viz., a yellow or purple cast, and had a lower contrast. 

The performance Index ranking averaged for each image product Is 
summarized by photographic parameters in Table l\,h. Here again, the current 
product was ranked the highest overall and the highest for all photo parameters 
except one, viz., saturation of the wheat class. 

A small number of anomalies occassional ly appeared on Product I and 
Hocutt,, which did not appear on the Kaneko and Kraus products. Inspection 
of the spectral data indicated that these anomalous pixels fell within the 
ranges defined for these two products and, therefore, were probably the 
result of system error. On the images that were inspected, these anomalous 
pixels wefe easily identified as such, so that in an operational system 
they should not confuse a ski 1 led analyst. The cause of these anomalous 
pixels, however, should be determined. 

Quantitative Evaluation 


The extent of pixel compression for each segment in the lower and 
upper band range is given in Tables A. 5a - A.5j and a summary of the average 
compression by segment is given in Table A. 6. Again, as had to be expected 
from the formulae, the most compression occurred In Product 1, followed 
by Kaneko and Hocutt. All products, however, minimized the compression 
In the lower range which is desired to maximize the detection of wheat 
early in its development. The spectral data for a small sample of wheat 
fields taken throughout the growing season from a number of spring and 
winter wheat segments indicated that none of these data were compressed 
in either the lower or upper ranges. Finally, it must be recognized that 
even in the extreme cases, the percentage of points that are compressed 
by any image type is very low, and little data are lost, regardless of 
the method used to calculate image specifications. 


Conclusions and Recommendations 
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Based on the qualitative and quantitative results of this limited 
study It has been concluded that the Product 1 Is the best currently 
available photographic product to be used by the A|'s. The Image quality 
of Product 1 was consistently good while that of the alternative products 
varied considerably, (Much of this variability may have been due to Improper 
film exposure or processing and not due to the manner In which the gains 
and biases were calculated. This could not be determined from the Infor- 
mation made available to RSRP personnel). Maintaining consistent, good 
Image quality Is extremely Important for an operational project such as 
LAC IE which Is dependent upon the processing of large volumes of A I data, 
so that analyst fatigue Is minimized. 

Of the three alternate products tested, Kraus Is probably the most 
logical one to use as an ancillary Image to Product 1, Although this 
conclusion Is not supported by the results of this study, It Is clear from 
the formulae used to calculate Input parameters to the PFC that the 
Kraus product Is the only one In which spectral compression of pixels Is 
Impossible. Therefore, on the Kraus product, all spectral data would be 
presented to the Al, thereby increasing the possibility pf early detection 
of wheat. 

It IS recommended that the development of ancillary Image products 
continue. This development, however, should have the objective of producing 
an image that emphasizes low density wheat only. It would not be necessary, 
therefore, to produce an Image that represents the full range of vegetation 
within the segment since this would be an ancillary product that would be 
used In conjunction with Product 1. This need for reduced spectral range 
should Increase the flexibility In the data manipulation to meet the 
desired objective. 


Table A*2 j 

The Segments and 
PFC Products 

Their Acquisition Dates Used to Evaluate Four 
• 

Segment 

Location 

Dates Available 

1178 

SE Kansas 

22 Nov 75, 10 Dec 75, 27 Mar 76, 3 May 76 , 20 May 76 

1 181 

SE Kansas 

11 Dec 75, 16 Jan 76 , 22 Feb 76 , 10 Mar 76, 



3 May 76 , 14 July 76 

1232 

SW Oklahoma 

31 Dec 75, 23 Feb 76 , 16 Apr 76, 22 May 76 , 24 June 76 

1521 

SW Minnesota 

23 May 76 , 10 June 76 , 21 Aug 76 

1538 

NE Montana 

25 Apr 76 , 13 May 76 , 31 May 76 , 18 June 76 , 5 July 76, 



23 July 76 , 11 Aug 76 

1574 

E. Central Nebraska 7 Nov 75, 30 Dec 75, 18 Jan 76, 4 May 76 , 9 June 76 

1622 

NE North Dakota 

19 Apr 76, 7 May 76 , 12 June 76 , 30 June 76 , 



19 July 76 , 6 Aug 76 , 24 Aug 76 , 10 Sept 76 * 

1681 

NE South Dakota 

6 May 76, 24 May 76, 10 June 76, 17 July 76 , 4 Aug 76 * 

1851 

NW Kansas 

19 J^n 76 , 6 Feb 76 , 24 Feb 76 , 13 Mar 76 , 31 Mar 76 , 



6 May 76 , 17 July 76 

1865 

SW Kansas 

15 Dec 75, 7 Feb 76 , 13 Mar 76 , 18 Apr 76 , 6 May 76 


No ground data available 
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Table A, 3: Ranking of PFC Products Based on Photographic Parameters 


SEGMEST 

PRODUCT 1 
PERFORMANCE 

RANK 

KANEKO 

PERFORMANCE 

RANK 

KRAUS 

PERFORMANCE 

RANK 

HOCUTT 

PERFORMANCE 

RANK 

1 181 

1.2 

1 

2.8 

3 

^f.O 

A 

2.1 

2 

1178 

1.8 

1 

3.0 

3 

3.1 

h 

2.1 

2 

1232 

1.3 

1 

3.1 

3 

3.5 

k 

2.2 

2 

1521 

2.1 

2.5 

1.7 

1 

2.1 

2,5 

4.0 

4 

1S3.3 

1.5 

1 

2.7 

3 

3.3 

4 

2.5 

2 

157V 

2.0 

1 

2.9 

3.5 

2,9 

3.5 

2.1 

2 

1851 

1.5 

1 

2.5 

3 

A.O 

4 

2.0 

2 

1861 

1.6 

1 

2.0 

2 

3.6 

4 

2.7 

3 

1865 

1.6 

1 

3.0 

3 

3.7 

4 

1 .8 

2 
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Table h.k: The Performance Index/Ranking of All Segments Summarized by the Photographic Parameters Used 

for Qualitative Evaluation of Four PFC Image Types 
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Table A. 5b: Pixel Compression for Segment ll8l — Lower Range (L) and Upper Range (U) 



Table =A. 5c: Pixel Compression for Segment 1232 — Lower Range (L) and Upper Range (U) 
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Table A.5d; Pixel Compression for Segment 1521 — Lower Range (L) and Upper Range (U) 
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Table A.5e ; Pixel Compression for Segment 1538 — Lower Range (1.) and Upper Range (U) 
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Table A.5f : Pixel Compression for Segment 157^ — Lower Range (L) and Upper Range (U) 
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Table A.5g : Pixel Compression for Segment 1622 — Lower Range (L) and Upper Range (U) 
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Table A.5h : Pixel Compression for Segment l68l — Lower Range (L) and Upper Range (U) 
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Table A.5i : Pixel Compression for Segment 1851 — Lower Range (L) and Upper Range (U) 
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Table A.5j ; Pixel Compression for Segment 1865 Lower Range (L) and Upper Range (U) 
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SEGMENT 


Table A. 6.; Summary of Average Pixel Compression — Lower Range (L) and Upper Range (U) for All Segments 
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Total Pixels Per Segment = 22,932 


